Mid-infrared plasmons in scaled graphene nanostructures 
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Plasmonics takes advantage of the collective response of electrons to electromagnetic 
waves, enabling dramatic scaling of optical devices beyond the diffraction limit. 
Here, we demonstrate the mid-infrared (4 to 15 urn) plasmons in deeply scaled 
graphene nanostructures down to 50 nm, more than 100 times smaller than the on- 
resonance light wavelength in free space. We reveal, for the first time, the crucial 
damping channels of graphene plasmons via its intrinsic optical phonons and 
scattering from the edges. A plasmon lifetime of 20 femto-seconds and smaller is 
observed, when damping through the emission of an optical phonon is allowed. 
Furthermore, the surface polar phonons in Si0 2 substrate underneath the graphene 
nanostructures lead to a significantly modified plasmon dispersion and damping, in 
contrast to a non-polar diamond-like-carbon (DLC) substrate. Much reduced 
damping is realized when the plasmon resonance frequencies are close to the polar 
phonon frequencies. Our study paves the way for applications of graphene in 
plasmonic waveguides, modulators and detectors in an unprecedentedly broad 
wavelength range from sub-terahertz to mid-infrared. 
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Introduction 

Plasmonic nanostructures allow the confinement of electromagnetic energy beyond the 
diffraction limit, enabling applications ranging from metamaterials, quantum optics, and 
photovoltaics to photodetectors and biological sensing 1 . In particular, graphene based 
plasmonics has garnered intense interests due to its unique tunability " and presumably 
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long plasmon lifetime ' and high degree of electromagnetic confinement ' . The 
electromagnetic response of patterned graphene micro-structures in the terahertz 
frequency regime was studied recently 4 ' 8 ' 9 , and was found to be dominated by the 
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plasmon excitations of Dirac fermions ' . Here we demonstrate that the dispersive nature 
of the 2-dimensional plasmon enables engineering of graphene plasmon resonances into 
the mid-IR regime by patterning graphene into nanostructures down to 50 nanometers. In 
contrast to their microstructure counterparts, whose plasmon resonances are in the 
terahertz frequency range, the mid-IR plasmonic response has been found to be affected 
strongly by interactions with substrate phonons and the graphene intrinsic optical 
phonons. These interactions lead to the renormalization of the plasmon dispersion and 
govern the lifetimes of plasmons. Moreover, an edge related effect on the mid-IR 
plasmonic responses is explicitly identified, which complements not only the extensive 
electronic transport studies performed previously on similar graphene nanostructures 10 " 12 , 
but also the research field concerning plasmons in metal cluster 13 ' 14 . 

Experiments and methodology 

Graphene nanostructures such as nanoribbon, antidot and nanodisk arrays are defined 
using electron-beam lithography on chemical-vapor-deposition (CVD) grown graphene 15 . 
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Different supporting substrates such as SiC>2 and DLC 16 are used on top of highly 
resistive silicon wafers. The mid-IR transmission measurement scheme is shown in Fig. 
la. The setup consists of an IR microscope coupled to a Fourier transform infrared 
spectrometer (FTIR), used in conjunction with a broadband polarizer. Details of the 
nanofabrication and infrared measurement are outlined in the Methods section. 

The qualitative behaviour of plasmonic response in graphene nanoribbons and nanodisks 
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are similar, except for a distinct polarization dependence in the former ' . In this paper, 
we focus on graphene nanoribbons, where a typical scanning electron micrograph is 
displayed in Fig. lb. The electromagnetic responses of these graphene nanoribbon arrays 
are characterized by their extinction spectra \-T per /To or l-T//To, where T per (T//) is the 
transmission through the ribbon array of the light with electric field perpendicular 
(parallel) to the ribbon, and T is the transmission through the substrate without graphene 
(but with SiC>2 or DLC). Fig. lc shows the extinction spectra of a ribbon array on SiC>2 
substrate, with widths of around 240 nanometers. The two curves correspond to incident 
light with polarization perpendicular (red) and parallel (grey) to the ribbons, respectively. 

Due to the excitation of localized plasmons, the extinction spectrum with perpendicular 
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polarization has prominent resonance peaks ' . The multiple resonance peaks observed 
here are in sharp contrast to the far-IR spectra of graphene microstructures, which usually 
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display a single strong resonance ' . Without plasmon excitations, the spectrum with 
parallel polarization lacks resonance peaks except for small features from 1000 cm" 1 to 
1200 cm" 1 . These features are related to the fast-varying dynamic dielectric function of 
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the SiC>2 substrate , and weak (< 1% in extinction) compared to the prominent plasmon 
peaks. In this mid-IR wavelength range, if the light polarization is parallel to the highly 
doped graphene ribbons (as reported in this work), their interaction is insignificant due to 
the Pauli blocking of interband transitions and the weakness of the free-carrier intraband 
transitions in such high frequencies . As a result, the transmission spectra of the ribbons 
for the parallel polarized light (7>/) are very close to the reference (To) and in subsequent 
presentations will serve as reference for the extinction spectra in the case of 
perpendicularly polarized light. In other words, the extinction spectrum is defined as 1- 
Tper/T// subsequently (See Methods). 

The multiple resonance structures observed in the spectra of Fig. lc are due to the 
interactions of graphene plasmons with the substrate phonons, as we will show in 
subsequent discussions. In fact, the Fuchs-Kliewer 19 surface-optical (SO) phonons are 
well-known in polar semiconductor surfaces and have been extensively studied in 
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conventional 2-dimensional electron systems ' . When graphene is placed on a polar 
substrate (e.g. SiCh) accommodating these SO phonons, the long-range Frohlich coupling 
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can mediate interactions with the electronic degrees of freedom in graphene ' , 
including the collective plasmon modes 24 . These interactions produce hybridized 
plasmon-phonon modes in the vicinity of their crossing energies and were recently 
observed by different techniques such as electron-energy loss spectroscopy in epitaxial 
graphene on SiC substrate and near-field optical nanoscopy of graphene on Si02 
substrate 5 . Before we present extended studies of these hybridized plasmon-phonon 
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modes, it is instructive to first consider a simpler system, i.e., graphene on a non-polar 
substrate DLC. 

Plasmons of graphene nanoribbons on DLC substrate 

Fig. 2a shows the extinction spectra of graphene ribbons with various widths Won a DLC 
substrate. Indeed, there is only one prominent plasmon resonance peak for each spectrum, 
unlike the spectra on SiC>2 substrate. Here, we would expect that the plasmon resonance 
frequency follows a simple yjq dispersion in the long wavelength limit, with q being the 
plasmon wave-vector, as predicted by the Random Phase Approximation (RPA) theory 
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for the linear response of a 2-dimensional electronic systems and graphene ' . For 
localized plasmons in graphene ribbons, q is simply given by n/W 29 . The inset of Fig. 2a 
shows the extracted plasmon frequencies from Fig. 2a as a function of n/W, where salient 

disagreement exists when compared to the simple scaling co pl ^Jn IW predicted by 

theory. This discrepancy cannot be due to the breakdown of the long wavelength 
approximation, because the wave-vector here is two orders of magnitude smaller than the 
Fermi wave-vector kf. Correction due to quantum confinement effect is also negligible for 
dimensions larger than 10 nm . 

Prior experiments with graphene nanoribbons have shown that it is often necessary to 
postulate an electrical width smaller than its physical width in order to explain their 
electrical transport behavior 10 " 12 . Indeed, fabrication of graphene ribbons by e-beam 
lithography and oxygen plasmon etching can impart atomic scale defects and complicated 
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edge chemistry that renders the edges electrically "inactive" . In a similar fashion, we 
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introduce an effective ribbon width defined as W e =W- Wo. Fig. 2b displays the plasmon 

frequency as a function of q = n/W e . The dashed lines are -^fittings, which now agree 

very well with the experiments for two different doping levels, if Wo = 28 nm is chosen. 
Using this methodology with the same Wo, we can also consistently explain the more 
complicated hybrid plasmon-phonon dispersions in graphene on SiC>2, to be discussed in 
the next section. 

Plasmons of graphene nanoribbons on Si0 2 substrate 

Fig. 3a shows the extinction spectra for graphene ribbons on SiC>2 substrate with W 
ranging from 60 to 240 nm, where we have labeled the three major resonance peaks and a 
weak higher order resonance 9 within our measured frequency range of 650 to 6000 cm" 1 . 
We identify four key features in the observed spectra: (1) Dispersion: all resonance 
peaks blue-shift as ^decreases but at very different rates. In particular, peak 3 disperses 
up in frequency at a much faster rate. (2) Spectral weights: The spectral weight transfers 
from peak 1 to peak 2 and 3 with decreasing W, and eventually peak 3 retains all of the 
spectral weight for W < 90 nm. (3) Linewidths: the resonance linewidth for peak 3 
increases with decreasing W, while those of peaks 1 and 2 remain almost constant. In fact, 
the resonance linewidth for ribbons on DLC discussed in previous section also exhibits 
the same trend as that of peak 3 in Fig. 3a. The linewidth is directly related to the 
plasmon damping and will be covered in detail in the following section. (4) Line shape: 
peak 3 is very asymmetric, especially for relatively narrow ribbons. As detailed in the 
Supplementary Information, the line shape can be well described by the Fano resonance 
model 33 . 
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The multiple resonances observed in the spectra are due to the plasmon-phonon coupled 
modes and can be described within a generalized RPA theory 24 ' 34 ' 35 . Directly relevant to 

our experiments is the loss function 3 [l/ £ r T pa ] , defined as the imaginary part of the 

inverse dielectric function, which describes the ability of the system to dissipate energy 
via plasmon excitations. Our calculations include both interactions with the relevant 
substrate SO phonons 17 at a> spl = 806cm 1 and co sp2 = 1168cm 1 and the graphene intrinsic 

optical phonon modes 36 at co = 1580cm 1 . Detailed description of the theory is provided 
in Supplementary Information. 

Fig. 3b shows an intensity plot of the loss function overlaid on the extracted resonance 
peak frequencies obtained from the spectra in Fig. 3a. We applied the same methodology 
used in previous section, i.e., defining the wave-vector of the localized plasmon to be q = 
k/(W-Wo), with the same Wo = 28 nm. Excellent agreement with the experimentally 
observed plasmon-phonon dispersions can be obtained using Frohlich coupling strength 
of F* pi ~ 0.2meV and F s 2 p2 ~ 2meV , where the latter includes the mode's double 
degeneracy. The Frohlich coupling quantifies the hybridization of the plasmon-phonon 
mode and their anti-crossing energy splitting and is defined as F?„ =^r( - ' - - ' )in 

the model (see Supplementary Information), where e high ( e low ) is the high (low) 

frequency dielectric constant of Si0 2 . The values extracted from fitting of the hybrid 
plasmon-phonon dispersions, although of the same order as quoted values used in the 
literatures ' , can serve as a better estimate. In addition, the second order dipolar 
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mode 9 ' 29 with wave-vector given by q = 3k/(W-Wq) also coincides with the dispersion of 
its fundamental mode. The above agreement further validates the premise that the 
effective width W-W defines the plasmon momentum, with a Wo = 28 nm consistently 

explains both plasmon and hybrid plasmon-phonon dispersions on both DLC and SiC>2 
substrates, respectively. 

Furthermore, we note that the loss function can also capture the plasmon peak intensity 
evolution for the three hybrid plasmon-phonon branches, featuring two anti-crossings and 
spectral weight transfer from the low frequency to high frequency plasmon branch (peak 
3) with increasing q. See also Supplementary Information. 

Origins of the Plasmon Damping 

Many plasmonic applications can benefit from a long plasmon lifetime 1 . However, 
depending on the system of interest, plasmon can decay into photons via radiative 
processes , into electron-hole pairs via inter- or intra-band Landau damping . In 
addition, inelastic scattering with phonons 14 ' 40 and elastic carrier scattering processes 9 ' 14 
can also contribute to the plasmon damping. The former two processes can be ruled out 
since contribution from radiative damping in graphene nanostructures is expected to be 

8 9 

negligible ' and that plasmon excitations in our devices lie outside of the Landau 
damping regime, as indicated in the Fig. 2b and Fig. 3b. However, the latter two 
processes can be important with large plasmon energies (i.e. %CQ pl >0.2eV) and small 

dimensions respectively. Lastly, interactions with the SO phonons can also influence the 
lifetimes of the hybrid plasmon-phonon modes, especially when the energy of the hybrid 
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mode is close to that of the SO phonons 41 . Our scaled graphene nanoribbons with 
different widths enable us to identify the damping mechanisms of graphene plasmons for 
the first time in a comprehensive manner. 

From the extinction spectra in Fig. 2a and 3a, the plasmon damping rate F p can be 
extracted from the resonance linewidth using a simple procedure as described in the 
Supplementary Information. Figure 4a plots the extracted r p as a function of the effective 
ribbon width W e . For the hybrid plasmon-phonon modes on Si02, only peak 3 is plotted 
since it exists for all of the different size ribbons and its damping rate evolves 
dramatically. We analyze first the extracted F p for DLC case. Since its plasmon 
resonance energy is below the optical phonon energy, the dominant broadening 
mechanism involves the elastic carrier scattering processes. Indeed, it's damping rate 
increases as W e decreases, indicative of scattering at the edges of the ribbon. Similar 
damping effect for plasmons in metallic nanoparticles has been extensively studied 13 ' 14 . 
This effect can be modeled using 13 ' 14 T p = T + a/W e , where F is the damping rate due 

to the carrier scattering in bulk graphene which was measured to be about 69 cm" 1 from 
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the Drude response of large area, unpatterned graphene nearby on the same DLC 
substrate and a is a fitting parameter. The same formalism was used to describe the 
damping in metal nanoparticles. Here we also want to emphasize that in micro-scale 
plasmonics structures, the damping rate r is almost identical to that extracted from 
Drude response, indicating the negligible impact of the edges in micro-scale structures. In 
the fit, a~2x\0 6 m/s , of the order of Fermi velocity, as expected from the edge 
scattering point of view. Lastly, we note that, in principle, the inhomogeneity of the 
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ribbon width may also lead to an apparent resonance broadening effect; however, 
estimates based on the SEM characterization, the doping dependence and ribbon width 
dependence of the plasmon peak broadening for the ribbons on DLC substrate, all suggest 
the negligible contribution of the inhomogeneity, see Supplementary Information. 

Despite similar characteristics of nanoribbons fabricated using CVD graphene (i.e. 
similar r and doping), the plasmon damping of ribbons on Si02 are different from its 
DLC counterpart as shown in Fig. 4a. First, we emphasize that the damping of hybrid 
plasmon-phonon modes also depends on its relative plasmon versus phonon character. 
For example, the mode described by peak 3 in Fig. 3 a is predominantly phonon-like for 
large W e since it resonates near the SO phonon frequency co sp 2. As a result, its damping is 
determined by the SO phonon lifetime instead (see Supplementary Information), which 
typically is in the picosecond regime. This explains the lower damping rate of this mode 
compared to the DLC case in the large W e limit, as shown in Fig. 4a for W e > 120 nm. 
Similarly, the small linewidths for peaks 1 and 2 in Fig. 3a can also be accounted for by 
the same reasoning. On the other hand, for narrower ribbons (i.e. W e < 90 nm), the mode 
described by peak 3 disperses rapidly in frequency and becomes predominantly plasmon- 
like. However, in contrast to the DLC counterpart, it has a plasmon frequency larger than 
the graphene optical phonon frequency of 1580 cm" 1 . Therefore, the plasmon can decay 
into electron-hole pair via the emission of a graphene optical phonon, which might 
explain the larger damping. This process is schematically illustrated in Fig. 4b and the 
following quantitative modeling verifies this mechanism. 
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Figure 4c plots the plasmon lifetime as a function of the plasmon frequency instead. The 
plasmon lifetime in graphene ribbons on DLC can be modeled by 



T ep can be calculated from the electron self-energy due to the interaction with the optical 

phonon °' 42 , albeit its effect is relatively small for our measured ribbons on DLC substrate. 
Lifetime modeling of the hybrid plasmon-phonon mode (i.e. peak 3) in the Si02 
counterpart is more involved, and we defer the details of the calculation to the 
Supplementary Information. After accounting for the plasmon decay channel via optical 
phonons, good agreement is achieved for plasmon lifetimes on both DLC and SiC>2, as 
shown in Fig. 4c. Here, we emphasize an important observation regarding the latter. In 
Si02, the hybrid mode begins with a long lifetime in the vicinity of the SP phonon energy 
at 1168 cm" 1 , but deteriorates rapidly to 20 femto seconds as its plasmon energy exceeds 
the optical phonon energy. This observation clearly highlights the important role of the 
various phonons on the plasmon lifetime. 

Lastly, we strengthen our case for the plasmon decay via optical phonons by studying the 
carrier density dependence of the plasmon linewidth. Figure 4d shows the extinction 
spectra of ribbons array with a width of 100 nm on SiC>2, at 4 different doping 
concentrations. The procedures for adjusting the doping concentration are described in 
the Methods section. We observed the increase of plasmon resonance frequency and 
amplitude with doping concentration, as expected from previously determined carrier 
density scaling law for plasmons in graphene . As the plasmon frequency exceeds the 
graphene optical phonon energy, we observed a substantial increase in its resonance 




where T is the additional damping due to optical phonons. 
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linewidth as presented in the inset of Fig. 4d. On the contrary, plasmon resonance 
linewidth decreases with increasing doping in the far-infrared regime for graphene micro- 
structures 8 . This observation further reinforces our picture of the plasmon decay channel 
via optical phonon emission. 

Outlook of the graphene plasmonics 

Through the scaling of plasmonic structures into the nanometer regime, we unravel the 
importance of both substrate and intrinsic optical phonons on the plasmon dispersion and 
lifetimes. These new findings have important implications to further our understanding of 
basic electronic properties such as carrier screening, scattering and energy 
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dissipation ' ' ' ' . Most importantly, our study provides important guidelines for 
future graphene -based plasmonic devices. Indeed, interaction with optical phonons in 
graphene limits the quality factor of the plasmon resonance beyond the optical phonon 
frequency. As a result, the optimal operation wavelength of future graphene plasmonic 
devices should lie between the sub-terahertz and 6.5 um. Nevertheless, the hybridized 
plasmon-phonon mode in the vicinity of the surface polar phonon frequency features a 
long lifetime. It might be promising to exploit this for high performance mid-IR photonic 
devices. 

Methods 

Large-scale graphene was grown using chemical vapor-phase deposition (CVD) on 
copper foil and then was transferred to Si02/Silicon or DLC/Silicon substrates 15 , before 
the subsequent e-beam lithography and oxygen plasma etching. The fabricated 
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nanoribbon arrays are designed to have width W equal to the spacing between the 
adjacent ribbons, corresponding to a filling factor of 50%. The array size is 70 |im by 70 
|im, significantly larger than the infrared beam size (-25 |im in diameter) used in the 
experiments. The as-prepared samples are hole-doped with a Fermi level of around -0.3 
eV. By exposing the samples to nitric acid vapor for 10 minutes, the Fermi level can be 
further decreased to below -0.5 eV. Baking the doped sample at 170 °C can reduce the 
doping significantly. As a result, by adjusting the doping and baking time durations, we 
can achieve various doping concentrations, which can be determined using the extinction 
spectra in the far-infrared . In the transmission measurements, we first recorded the 
transmission spectrum T per for the perpendicularly polarized light with respect to the 
ribbons, and the extinction spectra for this particular polarization {\-T per /To) is obtained 
using the transmission of the bare substrate (with oxide but no graphene ribbons) To as 
the reference. The extinction spectrum (l-T//To) for the light with polarization parallel to 
the nanoribbons was obtained in a similar manner. Light transmission through the 
substrate without graphene To has no polarization dependence. As shown in Fig. lc, the 
extinction spectrum for parallel light polarization (I-T//T0) was close to zero in the entire 
wavelength of interest. Consequently, we used \-T per /T// for the extinction spectrum of the 
light with perpendicular polarization for simplicity, since in this case no reference point is 
necessary. 
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Figure captions 

Figure 1. Schematics of the experiment 

(a) Mid-IR transmission measurement scheme of graphene nanoribbons. The excitation 
light is broadband. False colors are used for the excitation light (rainbow) and the on- 
resonance light in graphene nanoribbons (red), (b) A scanning electron micrograph of a 
typical array of graphene nanoribbons. The width of the ribbon is around 100 nanometers, 
(c) The extinction spectra of a ribbon array on S1O2 with the incident light polarization 
perpendicular (red) and parallel (grey) to the ribbons respectively. A weak 2 nd order mode 
is indicated. 

Figure 2. Plasmons in graphene nanoribbons on diamond-like-carbon (DLC) 

(a) The extinction spectra (\-T per /T//) of graphene ribbons on DLC with different ribbon 
widths. The spectra were referenced using transmission of light with parallel polarization 
(7}/). Inset: The plasmon resonance frequency as a function of the wave vector q = nlW, 
where ^is the width of the nanoribbon. Red curve is a fit according to Cfl p i ~y[q (b) The 

same data (red dots) as in (b), plotted as a function of wave vector q = 7t/W e , where W e = 
W-Wo is the effective ribbon width and Wo of 28 nm is used. The data from a lower Fermi 

level case is also plotted (grey dots). Dashed curves are fits based on the -yjq scaling. The 

shaded area is the intraband Landau damping region. The short bold black line just below 
1600 cm" indicates the graphene intrinsic optical phonon frequency. 

Figure 3. Plasmons in graphene nanoribbons on SiC>2 

(a) The extinction spectra (\-T per /T//) of graphene ribbons on SiC>2 with different ribbon 
widths. The spectra are vertically displaced for clarity. The vertical dash line indicates the 
graphene optical phonon frequency, (b) Plasmon frequency as a function of the wave- 
vector q (= 7vW e ) for peaks 1, 2, 3, and a 2 nd order dipolar mode peak in (a). The 
calculated loss function is plotted as a 2-dimensional pseudo-color background. The 
dashed line represents the plasmon frequency without considering the plasmon-phonon 
hybridization. Two surface polar phonons and the intrinsic optical phonon frequencies are 
indicated. See Supplementary Information for details of modeling. 

Figure 4. Origins of the plasmon damping 

(a) The damping rates of the plasmons of graphene ribbons on DLC (red) and SiC>2 (grey) 
as a function of the effective width W e . The data points are obtained from the linewidths 
of the spectra in Fig. 2a and Fig. 3a. For ribbons on SiC>2, peak 3 is used. The red dashed 
curve is a fitted curve based on the model described in the main text and the grey dashed 
curve is a guide to the eye using r °c (W -W o y x . The data points for ribbons on SiCh 

with resonance frequency above the optical phonon frequency are indicated by the yellow 
rectangle, (b) An illustration of the plasmon damping process through the emission of an 
optical phonon together with the creation of an electron-hole pair. Interband and 
intraband damping regions are indicated, (c) The plasmon lifetimes of ribbons on DLC 
(red dots) and S1O2 (grey dots) as the functions of the plasmon resonance frequencies. 
Dashed curves are calculated results, with the orange curve representing the contribution 
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from plasmon damping via phonon emission, the blue curve for contributions from both 
phonon emission and bulk scattering, and the grey (for SiCh) and red (for DLC) curves 
for contributions from phonon emission, bulk scattering and edge scattering. See 
Supplementary Information for details of calculation. The yellow region represents the 
frequency range larger than the graphene optical phonon frequency, (d) The extinction 
spectra of a ribbon with W of 100 nm on SiCh at four different doping levels. The inset 
shows the extracted full width at half maximum (FWHM) as a function of the 
corresponding plasmon resonance frequency. 
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Figure 1 Schematics of the experiment 
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Figure 2 Graphene nanoribbons on diamond-like-carbon (DLC) 
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Figure 3 Graphene nanoribbons on Silicon Dioxide (Si02) 
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Figure 4 Origins of the plasmon damping 
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I. CALCULATING THE RPA DIELECTRIC FUNCTION AND LOSS FUNCTION 



This section describes the modeling of the loss function in graphene on SiC>2 as shown in 
the intensity plot of Fig. 3b in the main manuscript. 

A. Intrinsic phonons 

The only intrinsic phonons with momenta and energies similar to the graphene plasmons 
in our experiment are the longitudinal/transverse optical (LO/TO) phonons near the T 
point, with energies hu op ~ 0.2 eV. On symmetry grounds jT], their coupling to electrons can 
be written as [21 [3] 

v f\ I u y (r) + iu x (r) \ j,, 

^e-op(r) = g \ = — j(r) x u(r) (1) 

y Uy (v) - iu x (r) J V F 

where vf is the Fermi velocity and go the coupling constant. The coupling constant can be 
estimated from the change with bond length of the hopping between nearest neighbor carbon 
7r orbitalsjU |5], go ~ dt/dl. The electronic Hamiltonian is described within each valley 
(and spin) in terms of the amplitudes on A/B sublattices, j(r) is the single-particle current 
operator and u(r) is the relative displacement of the two sublattices. Their representation 
in terms of electron and phonon ladder operators, i.e. and 6 q respectively, are given by, 

j>) = \ E W*' = i Ei q e iqr (2) 

q 




where a are the Pauli spin matrices and e q A are the polarization vectors. A denotes the 
phonon modes, p m is the mass density of graphene and A is its area. Using standard pertur- 
bation techniques, the effective electron-electron interaction mediated by optical phonons 
can be written as, 

Y&-* = ^ E ^ \ M ° P \ 2 ■ U = ^2 E V,Ajq • j-q (4) 

qA F qA 

where the scattering matrix elements and the free phonon Green's function are 

' M -l 2 = ^- ' V °^ = ,(( — ^ ( 5 ) 

2p m w op h((u + th/T op y - - 1 

where r op phenomenologically describes the phonon lifetime. 
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B. Surface polar phonons 



Polar substrates, such as SiC>2 and BN have optical piezoelectric modes at energies huj sp . 
These modes induce electric fields which couple to the carriers in graphene[Hl U\- At long 
wavelengths, the effect of these can be described in terms of the dielectric function of the 
substrate, 

Ke-sp =jYl M ^4+ q a k(& q A + bl qX ) (6) 

kq 

with the scattering matrix elements defined as, 

\M SV \ 2 = T , T 2 = —^( (7) 

^env ^low * ^env J 

where z is the graphene-substrate separation, T 2 describes the Frohlich coupling strength, 
tiow (thigh) are the low (high) frequency dielectric constant of the dielectric and e env is that 
of the environment. The effective electron-electron interaction mediated by surface optical 
phonons calculated from standard perturbation techniques yields [8], 

V !el = E l M ^! 2 V »P^ = E ^,APqP-q ( 8 ) 
qA qA 

where p q = Ek^k^k+q anc ^ ^xi 00 ) contains also a phenomenological phonon lifetime of r sp . 



C. Dielectric response 

The plasmon response of graphene begins with finding the dielectric function. A satisfac- 
tory approximation can be obtained by adding the separate contributions independently. An 
effective interaction between electrons is given by the sum of the direct Coulomb interaction 
v c{q) = e 2 /2geo and the two electrons interaction mediated by surface phonon v sPt \(q,uj). 
The RPA expansion of the dielectric function, e^ a (q, u), can be expressed with this effective 
interaction [91 [10] 

v eff {q,UJ) - -r - (9 J 



e r ^ a (q, u) 1 - [v c (q) + *£x v*p,a]II° (q, 



where U° pp (q,u) is the non-interacting part (i.e. the pair bubble diagram) of the charge- 
charge correlation function given by a modified Lindhard function [TT1 IT2] . 

<,(«>») = -T^Y, [ F w (k,q) (10) 

P,P (27r) 2 JBZ £k - £k+q + fiW + 2ft/ r e 



28 



where rap(£k) is the Fermi-Dirac distribution function, F nn >(k, q) is the band overlap function 
of Dirac spectrum and r e is the lifetime of electrons. While the polar surface phonons couple 
to the charge density operator, the intrinsic optical phonon couple instead to the current 
operator. Its contribution to the dielectric function is given by v op (q,co)H^Aq,cj), where 
Hjj(q,oj) is the current-current correlation function. We note that from the usual charge 
continuity equation, idtpq = q • j q , it follows that, 



qU j j(q,u)=uU Ptft {q,u)-v F { q-j q ,p_ q ) (11) 



where the second term in Eq. 11 is purely real and oc q 2 as calculated in Ref. [13]. The 
imaginary part of Ujj(q,cu) can be obtained just from x U PtP (q,cu)}. Collective modes 
with self consistent oscillations of the carrier charge can be obtained from the zeros of the 
full dielectric function 

e^ a (q, u) = t env - v c U° p p (q, u) - e env ^ v s P ,x^° p , p (q, w) - t env v op R° j:j (q, u) (12) 

A 

where e env is the dielectric constant of graphene's environment. 



D. Loss function 



Our spectroscopy experiments measure the extinction spectra defined as Z = —5T/Tq 
with 5T = T — T , where T (T ) is the measured transmission with (without) plasmon 
excitations. In the experiment, a superlattice of graphene ribbons of width W defines 
the momentum i.e. q = ir/(W — Wo). Wq accounts for the difference between physical 
and electrical device's width. Varying the frequency of the incident light excitation, u, 
polarized perpendicularly (and parallely) to the ribbon, allows one to quantify the extinction 
spectra Z(q,u), as first demonstrated in Ref. [13] . Resonance peaks in Z(q,u) corresponds 
to enhanced optical absorption by graphene originating from plasmon oscillations [T3H16] 
and can best be described by, 

(13) 

where the latter is known as the loss function, which describes the ability of the system to 
dissipate energy via plasmon excitations and can be calculated from Eq.[T2j 



Z(q,u) ~ — S 



rpa 
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Using the above theory, we plot the loss function in graphene on SiC>2 as shown in 
Fig. 3b of the main manuscript. The calculations include interactions with the intrinsic 
and Si0 2 substrate phonons. Graphene doping is assumed to be Ej = — 0.43eV and an 
effective e env = 1.5. The frequencies of the various phonon modes are assumed to be at 
uj op = 1580 cm -1 , u sp i = 806 cm -1 and u sp 2 = 1168 cm -1 . The damping time used in 
those plots are r e = 0.1 ps, r op = 70 fs and r sp = lps. The coupling parameters used are 
go = 7.7eVA -1 , x = 0.2 meV and J-^ p2 = 2meV. Note that another substrate phonon at 
u op = 460 cm -1 was not included in the calculation, given that our experiment data are far 
above that frequency. 

II. CALCULATING LIFETIMES OF THE PLASMON AND COUPLED 
PLASMON-PHONON MODES 

This section describes the modeling of the plasmon lifetime in graphene on DLC and 
Si02 as shown in Fig. 4a and 4c of the main manuscript. 

In the above previous analysis, the damping mechanisms for the plasmons are not dis- 
cussed. Exchange of energy and momentum during scattering of plasmons can bring it into 
the Landau damping regime, leading to finite damping. In fact, when the plasmon energy 
exceeds the optical phonon energy, it can decay into a phonon together with an electron-hole 
pair, in such a way that the total momentum is conserved. In a phenomenological way, 
this decay can be accounted for through the single particle excitations, which have a finite 
lifetime r e , when their energies exceed the optical phonon energy [T7], for example. Damping 
related to scattering with the ribbon's edges and a background damping due to impurities 
in the bulk can also be incorporated in r e . Finite phonon lifetime, r sp , can also influence to 
plasmon damping in the coupled plasmon-phonon modes. Below, we present our description 
of plasmon damping in the presence and absence of coupling with the surface phonon modes. 

We are interested in the regime where u> > vpq and Ej 3> hu. In this limit, 



where 5 e is the single-particle related damping in graphene defined as 5 e = l/r e . In the 




(14) 
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absence of substrate phonon interactions, such as the case of graphene on a DLC substrate, 
the plasmon frequency is simply u = u p i + i6 e , where u pl = q\E f\e 2 /2nh 2 e e env and 5 e also 
corresponds to plasmon damping. Note that in the regime we are considering i.e. u > vpq 
and Ef ^> fad, Landau damping is excluded. Guided by experiments, the plasmon hybridizes 
strongly with one of the surface phonon modes with u sp ~ 0.145eV on SiC>2 substrate. For 
u> > u sp , we can write a simpler dielectric function, 



rpa ^ 



(u + i8 e ) 2 (u + i5 sp ) 2 - id 2 p + &%, 



~ 2 i '4tt _ 

—UspJ 72 (15) 



where 5 sp is the surface phonons damping rate defined as 5 sp = l/r sp . The frequencies of 
the coupled plasmon-phonon modes can be obtained by setting e]! a = i.e., 

u 4 + i2u 3 (S sp + S e ) - uj 2 (lo 2 p + u p i{q) 2 ) - i2u(u 2 p S e + uo pl {q) 2 5 sp ) + u pl (q) 2 (u 2 p - u 2 p ) = 

(16) 

which can be solved numerically. In the limit where 5 e = S sp = 0, it reduces to a simple 
biquadratic equation with coupled plasmon-phonon modes solutions given by, 



U ± = —2— ± 2 (17) 

In the general case where 5 e = 5 p h ^ 0, we solve for the coupled plasmon-phonon modes 
via Eq. 16 numerically. However, in the q = limit, it can be shown by setting = in 
Eq. 15 that u = u sp — i5 sp . Therefore, the lifetime of the plasmon with frequency in the 



vicinity of the surface phonon frequency is determined by the surface phonon lifetime instead. 



In this work, we assume that r sp is constant, to be fitted to experiment. Here, we discuss 
model description of the electron lifetime r e . Including relevant scattering mechanisms in 
our experiments, r e is given by, 

T e (q, u) w [tq 1 + T edge (qY l + T ep (uoY l ] ^ (18) 

where r describes a background damping due to scattering with impurities and r edge (q) w 
ajiyV — Wo) b is related to scattering off the ribbon edges, tq ~ 85 fs as measured from the 
Drude response of large area, unpatterned graphene. a m 2 x 10 6 , of the order of Fermi 
velocity and b = 1 as discussed in the main text. r ep (u) is electron lifetime due to scattering 
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with optical phonons. It is related to the electron self-energy S ep via r ep = h/2^s[T lep ] given 



E ep (cu) 



-k B T \M op \ 2 V x {u x )°g°{k f + q,u + u x ) 



(19) 



where Q° is the electron Green function and ky is the Fermi wavevector. According to 
density functional calculations, the imaginary part of S ep can be approximated by[18j, 

1 " 



3[£ ep (u;)] = 70 \ fajJ + hw + Ef \ x 



erf| fa :^ Uerf r:^ l+2 



(20) 



where 7 describes the effective e-ph coupling and A ph accounts for various energy broad- 
ening effects such as the deviation from the Einstein phonon dispersion model. They are 
estimated to be 70 ~ 0.018 and A p h ~ 50meV from density function calculations |18j. 



As discussed previously, in the absence of interaction with the surface phonons, the 
plasmon lifetime is simply r e (q,u = oj p i). In the presence of interaction with surface 
phonons, the plasmon lifetime for the plasmon-phonon coupled modes can be solved via 



Eq. 16 numerically, with 5 e ~ [r e (q, u = u) p i{q))]~ 1 . The computed plasmon damping 
rates or lifetimes on DLC and SiC>2 substrates are shown in Fig. 4a and 4c of the main 
manuscript. There, we assumed graphene doped at Ef ~ — 0.43eV. 



Previously, we also show in Fig. 3b of the main manuscript the RPA loss function in SiCV 
However, the electron lifetime there was simply assumed to be constant. In Fig. S[TJ we 
calculate again the RPA loss function in SiC>2, but this time including the electron damping 
5 e which describes our experiment as detailed in this section. After the inclusion of a more 
accurate description of electron damping, we note that the loss function can capture very 
well the plasmon peak intensity evolution for the three hybrid plasmon-phonon branches as 
seen experimentally, featuring two anti-crossings and spectral weight transfer from the low 
frequency to high frequency plasmon branch (peak 3) with increasing q. 



III. RESONANCE LINESHAPE AND THE EXTRACTION OF THE PLASMON 
DAMPING 

This section describes the extraction of the plasmon damping rate from the measured 
extinction spectra, used for Fig. 4 of the main manuscript. 
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The far-IR plasmon resonace lineshape of graphene micro-disks and ribbons can be well 
described by a damped oscillator modeipiJ], which is derived from Drude conductivity. In 
the mid-IR regime, the lineshape becomes asymmetric, as seen for the third peak in Fig. 3a 
of the main manuscript. Because the lineshape of the resonance peaks in the vicinity of the 
substrate phonon frequencies might be affected by the plasmon-phonon hybridization, here 
we focus on the third peak of relatively narrow ribbons on SiC>2 substrate, whose resonance 
frequency is far away from those of substrate phonons. Fig. ^2] shows a spectrum (black 
curve) of a ribbon array with width W = 85 nm. The spectrum is very asymmetric and a 
Fano resonance model[20l I2T] can well describe it, as shown by the red curve in the figure. 
In the Fano framework, the extinction spectrum is expressed as 

T per = 2p (q f + T]) 2 

T„ nT p (l + q 2 ) (1+ V 2 ) 1 ] 

where p is a parameter for the amplitude, qf is the Fano parameter, T p is the plasmon 
damping rate, and rj = 2(ui — uo)/T p with uq being the center frequency. This equation is 
used to fit the spectra shown in the main manuscript to extract the plasmon damping rate 
T p . Typical values for Fano parameter qf are around 3. 

This Fano type resonance indicates that the plasmon resonance is interfered by a broad 
background continuum. As demonstrated before, the optical conductivity of graphene in 
the far-IR is Drude like. However, in the mid-IR range we are dealing with here, it's in the 
Pauli blocking regime, i.e., the optical conductivity has a very weak Drude response tail and 
some residue conductivity due to many-body effect. It is this residue conductivity serving 
as the broad continuum to form a Fano type resonance with the plasmon excitation. 



IV. RESONANCE BROADENING EFFECT DUE TO THE RIBBON WIDTH 
INHOMOGENEITY 

Long wavelength variations in the ribbon's width can also lead to an apparent reso- 
nance broadening effect indistinguishable from lifetime broadening effects. In DLC rib- 
bons, the plasmon dispersion is described by a simple u pl = aq dispersion, where a = 
\Ef\e 2 /2irh 2 eQe env . In this case, the resonance broadening 6u p i associated with a character- 
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istic width inhomogeneity of A is given by, 

5oo pl « -^-^ (22) 
2Wi 

In comparison to the plasmon lifetime broadening due to scattering off the edges, it also has 
a power law behavior of the form a/W^, the scaling exponent b in this case is | instead of 
1. Fig. £[3] plots the experimentally extracted plasmon damping as a function of the effective 
width W e for two different doping levels. Least-square-error fit to the data yields b = 1.0, 
indicating that the broadening is due to carrier scattering off the edges. Furthermore, 
as shown in Fig. £[3j there is also no noticeable dependence on doping, where the width 
inhomogeneity effect would suggest otherwise. 
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FIG. S 1: Graphene loss function on Si02 substrate. Calculated RPA loss function 
9? [l/ e^? a ] , including interactions with the intrinsic and Si02 substrate phonons. Graphene doping 
is assumed to be Ef = — 0.43eV and an effective e env = 1.5. Shaded regions represent the intra- 
band Landau damping regime i.e. huj/Ef < q/kp. Dashed line on the left plot is calculated from 
the classical plasmon dispersion = e 2 qvFkF/(2Trheoe env ). The frequencies of the various phonon 
modes are assumed to be at co op = 1580 cm -1 , lo sp i = 806 cm -1 and u sp 2 = 1168 cm -1 . The lifetime 
associated with the phonons used in these plots are r op = 70 fs, T sp i = 0.5 ps and T sp \ = 0.2 ps. 
Calculations include damping of single particle excitations 6 e as described in the Suppl. info. text. 
The coupling parameters used are go = 7.7eVA -1 , F^ pl = 0.2 meV and Tl p2 = 2meV. 
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FIG. S 2: Resonance lineshape and fitting. Measured extinction spectra of a W = 85 nm 
ribbon array on SiC>2 substrate, where we show only part of the spectra relating to the peak 3 as 
described in the main manuscript. The data can be fitted well by a Fano model, with model and 
parameters described in Suppl. Info text. 
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FIG. S 3: Damping rate scaling with width. Measured plasmon damping rate in graphene 
ribbons on DLC substrate as a function of the effective width W e = W — Wo, where W is the 
physical width and Wo = 28 nm. See also main text. The plasmon damping rate is best described 
by the scaling relation To + a/W^, where To = 69cm -1 is related to the background damping, a 
quantity determined from the Drude response of large area graphene. a and b are obtained from 
least-square-error fit to the experimental data, b was found to be 1.0. 
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